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Abstract 

For triples of probability measures, pure quantum states and mixed 
quantum states we obtain the exact constraints on the fidelities of pairs 
in the sequence. We show that it is impossible to decide between a 
quantum model, either pure or mixed, and a classical model on the 
basis of the fidelities alone. Next, we introduce a quantum three state 
invariant called phase and show that any sequence of pure quantum 
states is uniquely reconstructible given the fidelities and phases. 



1 Preliminary 

The transition probability or overlap between states is one of the fundamental 
ingredients of quantum mechanics. It lies at the basis of its probabilistic 
interpretation. Two states with an overlap close to one are almost equal. 
This motivates the term fidelity which will be used as a synonym of transition 
probability. 
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Also in the context of classical probability there is a quite common notion 
of distance between probability measures, the Hellinger distance based on 
overlaps of densities. This distance is given by djj(/x; A) = 1 — A(fi; A), where 
A is called the affinity between A and \x. Affinity is the classical counterpart 
of transition probability. 

We shall be concerned here with the notion of fidelity in three different set- 
tings: probability measures, pure quantum states and mixed quantum states. 
As there are natural inclusions between these three sets extending the affinity 
between classical measures to fidelity for mixed quantum states, we gradually 
widen our scope when passing from probability measures to mixed quantum 
states. 

Notions of closedness of states belonging to a sequence generated by a dy- 
namics or by a coding procedure are obviously important to decide either on 
the regular or the chaotic nature of the sequence [3] or on compressibility 
issues, see [7, 6] . In previous work, we used the spectrum of the Gram ma- 
trix of the sequence for this purpose, see [2, 4]. This spectrum is however a 
complicated function depending on the full sequence. 

The point of this paper is twofold. In Section 2, we analyse sequences of 
three states. This is the simplest non-trivial situation. We obtain the precise 
constraints on the fidelities of pairs of states in the sequence for measures 
and for pure and mixed quantum states. They turn out to coincide. In other 
words, on the basis of an admissible triple of fidelities it is impossible to 
decide whether one deals with probability measures or with pure or mixed 
quantum states. 

Next, we introduce in Section 3 a quantum three state invariant called phase 
and we show that, for pure quantum states, any arbitrary sequence can be 
uniquely reconstructed on the basis of the fidelities of pairs and phases of 
triples in the sequence. 

We recall the basic definitions. 

The fidelity between two probability measures A = (Ai, A2, • • • , A^) and /x = 
/i2, ■ ■ ■ , A*jv) on an event space with N elements is given by 
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We may extend fidelity to general probability measures by using the Radon- 
Nikodym derivative. 
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For quantum systems, we call expectation functionals on the observables 
states. E.g., in standard quantum mechanics a normalised vector ip in the 
Hilbert space S) of the system yields a pure state 

P<p{X) := (<f,X if), X an Hermitian operator on $). 

Writing 

P v {X)=Ti(\<p)(<p\X) 

we see that P v is uniquely determined by the one-dimensional projector on 
the subspace Ctp and we shall therefore identify |y?)(y?| with P^. A (normal) 
mixed state on S) is then a general density matrix p. 

The fidelity between P v and P^ is the usual quantum mechanical transition 
probability between ip and ip 

F(P V -P^) := \(<p,ip)\ 2 = Tr (P v P+) . (2) 

For mixed states, Uhlmann extended the notion by using the purification 
procedure, i.e. by obtaining a mixed states as a marginal of a pure state on 
a composite system, see [8]. This leads to 

F(p ; a) := (Tr \J ' p^ a p^j = ^Tr \J~o^~pai^ . (3) 

An isometric transformation U induces the transformation p i— > U pU* = 
poAd(U*) on density matrices. The spectrum of p, taking degeneracies into 
account, remains invariant up to multiplicities of zero. This is also the case 
for the fidelity 

F(p;a) = F(U pU* ;UaU*) = F(p o Ad(CT) ; a o Ad(ET)). 

Again, fidelity between states can be extended to general quantum proba- 
bility spaces, see [1, 8]. It is straightforward to check that the three defini- 
tions (l)-(3) of fidelity are compatible. 

A simple example is given by states on the matrices of dimension 2. In this 
case there is an affine isomorphism between the density matrices of dimension 
2 and the unit ball in R 3 explicitly given by the Bloch transformation 

p x := i(l + x ■ tr), x ER 3 , \\x\\ <1. 

Here cr are the usual Pauli matrices. The two dimensional case is rather mis- 
leading as the topological and convex boundaries of the state space coincide, 
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they are namely the unit sphere in R 3 . For the ci-dimensional matrices, the 
state space has d 2 — 1 real dimensions, its topological boundary d 2 — 2 while 
the pure state space has dimension 2(d — 1). So, the topological boundary 
has many flat pieces. The Uhlmann fidelity can be readily computed for the 
two dimensional case and one obtains 

F(p x ; Py )= l -(l + x-y + Vl-IWlVl-IMI 2 ) • 

In particular, for pure states, i.e. ||£c|| = \\y\\ = 1, the fidelity is the square 
of the cosine of half the angular distance between x and y. 



2 Fidelities of triples 

In this section, we consider a triple of states or measures and compute the 

fidelities pairwise. We then look for the constraints between the three fideli- 

iii 

ties F 12 , F 13 and F 23 . Plotting (F* 2 , F^ 3 , F 2 \) in R 3 we obtain subsets of (R + ) 3 
when the states run either through the measures, the pure or the quantum 
mixed states. These sets are obviously ordered by inclusion. Moreover they 
are closed, convex and compact. The convexity follows from the direct sum 
construction. E.g. for the Uhlmann fidelity 




= \ (Tr V Pi ^ Pi + Tr V Pi ^ pl) ■ 

Lemma 1. Let (p 1} <p 2 and ip 3 be three normalised vectors in a Hilbert space 
$) and let F jk := F(P LPj ;P ipk ), then 

F12 + Fi 3 + F 23 < 1 + 2^F 12 F 13 F 23 . (4) 

Proof. The three dimensional matrix 

/ 1 (y?i,y? 2 > {<Pi,<P3)\ 

\(<P3,<Pl) (V^3,^2) 1 / 

is positive semi-definite. In particular, det(G) > or 

< l + 2Ke(jy?i,y?2)(y?2,y?3)(y?3,y ? i)) - I (v?i,v?2)| 2 - I (v?i, V^3> I 2 - 1(^2, ^s)l 2 - 
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Hence 

F 12 + F 13 + F 23 = |(v?i,v? 2 >| 2 + |(y?i,y? 3 >| 2 + |(^2,y?3>| 2 

< 1 + 23?e((</?i,</? 2 ) (f2,f3) (^3,^1)) 

< l + 2\(<Pl,<P2)\\(<P2,<P3)\\(<P3,<Pl)\ 

= 1 + 2 a/ -F12-F13F23, 
which proves (4). □ 

Let us denote by C 3 the compact convex subset of R 3 

{x = (xi, x 2 , X3) I < Xj < 1 for j = 1, 2, 3 and x x +X2 +X3 < 1 + 2x1X2X3}. 

The convex boundary of C3 is 

C (C 3 ) ={(0, 0, 0), (0, 0, 1), (0, 1, 0), (1, 0, 0), (1, 1, 1)} 

U {x I < Xj < 1 for j = 1, 2, 3 and x\ + x\ + x\ — 1 + 2xiX 2 x 3 }. 

ill 

Lemma 2. C3 = {(^12,-^13,-^23)} where the fidelities are computed for prob- 
ability measures on a space with six points. 

Proof. The extreme point (0, 0, 0) can be realised with three degenerate mea- 
sures on a space of three points, (1,0,0), (0,1,0) and (0,0,1) on one with 
two points and (1, 1, 1) with a single point. Using convexity and the shape 
of C3, it is then sufficient to show that any extreme point of the surface 
{x I < xj < 1 for j = 1, 2, 3 and x\ + x\ + x\ = 1 + 2x1X2X3} can be 
obtained by transition probabilities between measures on a space with two 
points. 

Suppose that x x < x 2 < x 3 and put Ai := (1,0), A 2 := (xf, 1 — x\) and 
A 3 := (x 2 , 1 — x 2 ). We then have that 

F(\ 1 ;\ 2 )=x 2 1 , F(Ai;A 3 )=x2, and 



F(A 2 ; A 3 ) = (x x x 2 + y / (l-x?)(l-x2)) 2 . 

Solving the equation x\+x\ + x\ = I + 2X1X2X3 for X3 and taking into account 
that < Xi < x 2 < x 3 < 1 we find that F(\ 2 ; A 3 ) = x\. □ 
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Lemma 2 shows that any triple of transition probabilities that can be ob- 
tained by transition probabilities between pure quantum states can also be 
reached by classical measures. Our next aim is to show the rather surprising 
fact that we still remain in the set C 3 when we consider Uhlmann fidelities 
between three arbitrary mixed quantum states. 

Proposition 1. Let p±, p 2 and p 3 be three density matrices on a Hilbert space 
ft and let Fjk : = F(pj ; pk), then 

F\2 + -F13 + -F23 < 1 + 2 a/ -Fl2-Fl3-^23- 



Proof. We may assume that the pj have full support. The general case is 
obtained by continuity. Denoting F(px;p 2 ) 5 by a, we certainly have that 
< a < 1. Therefore 

< F( Pl ; p 3 ) + F(p 2 ; p 3 ) - 2a v / F(p 1 ; p 3 ) F(p 2 ; p 3 ) 

and we wish to obtain the upper bound 1 — a 2 . Using the polar decomposition 
and the assumption on the support of the pj, there exist unitaries U and V 
such that 

Pi Pi =U\p{ Pl I and p 2 2 pf = V |p 2 2 pf | . 

In these formulas \X\ denotes the absolute value of an operator X, i.e. \X\ := 
(X*X)2. We now express the fidelities as follows 



F 3 ( Pl ; p 3 ) = Tr V p\ Pi p\ = Tr|pf pf | = Tr U* pj pf . 

Using the Hilbert-Schmidt scalar product between Hilbert-Schmidt opera- 
tors, (X, Y)us '■= TrX*Y, the fidelities become 

(pi ; p 3 ) = (/, h)ns and (p 2 ; p 3 ) = h) HS 

with i i i 

/ := Pl U, g := p 2 2 V and /i := p| . 

We can then verify the following properties 

(f,h) H s = \(f,h) H s\, (g,h) H s = \(g,h)j3s\, 
|hs = II^IIhs = II^IIhs = 1, and 



\(f,g)\ = \TtU*p!pIV\ < sup \Trplp*W\ = F*(p l]P2 )=a. 

W unitary 

The statement of the proposition now follows from Lemma 3. □ 
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Lemma 3. Let f and g be normalised vectors in a Hilbert space ft and let a 
be such that \ (f, g)\ < a < 1, then 



sup 

h, \\h\\<l 



\(f,h)\ 2 +\(9,h)\ 2 -2a\(f,h)\\(g,h)\) = (l-a)(l+|</, <?») < (1-a 2 ). 



Proof. As the supremum is always non-negative, we may impose the addi- 
tional restriction h G span({/, g}). If h does not belong to this subspace, 
decompose it into h\®h 2 with hi G span({/, g}). Next replace h by /ti/||^i||. 
Evaluating the functional with this new h will return a value at least as large 
as that with the original h. 

So let h = af + (3g with a, (3 G C such that 

\\hr = \a\ 2 + W + 2Kz{aP(f,g)) = l. (5) 

Using this normalisation condition we compute 

K/,^)! 2 = |« + /3</,^>| 2 = i - |/?| 2 (i - |</,^>| 2 ) 

and 

\{g,h)\ 2 = \a{f,g) +]3\ 2 = 1 - |a| 2 (l - |(/, g)\ 2 ). 

Hence, the functional of h we have to maximise does not depend on the phase 
of a(3{f, g). The normalisation condition (5) can be satisfied if and only if 

M 2 + |/?| 2 -l| <2|a||/3||</,<?>|. 
Putting A := \a\, \i := and t := \ (f,g) \ we have to compute 

/ := sup (2 - (A 2 + /i 2 )(l - t 2 ) - 2ay/l - A 2 (l - t 2 ) y/l - /i 2 (l - t 2 )) 

subject to the constraints 

< A, 0</i, and |A 2 + fi 2 - 1| < 2Xfd 

with t satisfying < t < a < 1. The supremum is attained choosing A = /i, 
with A such that 

1 1 

< A < 



2(1 + t) ~ ~ 2(1 -*)" 
We obtain for I the value 

J = 2(l-a)(l-i T ^) = (l-a)(l + t)<(l-a 2 ). 

□ 
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3 Phase, a three state invariant 



For three pure states Uj := P (fj with tpj G and \\(pj\\ = 1, j = 1,2,3 we 
introduce the complex number e **( w i;^2;w 3 ) by the relation 

(<Pi, <f2) (<P2, <P3) (<P3, <Pi) = e l * (wi;W2;W3) (v?i, <p 2 ) (<p 2 , tfis) (<P3, <Pi 

= e^^F^u, ; cu 2 )F^(cu 2 ; u 3 )F*(u 3 ; cu 1 ). 

This definition implicitly assumes that all fidelities are strictly positive. It is 
easy to see that there is no continuous extension to the case of zero fidelities. 
Obviously $ is invariant under isometric transformations and 

$(cji; u 2 ; uj 3 ) = $(^ 2 ; ^3; ^1) and ®(u)i, uj 2 ; u 3 ) = -$(cj 2 ; Ui] uj 3 ). 

The relevance of other invariants than the fidelity was e.g. underlying the 
study of the spectrum of the Gram matrix of a sequence of states as in [2] . 
An other instance was considered in [5] where compressibility was shown to 
depend monotonically on the elementary symmetric invariants of the Gram 
matrix. Here we shall argue the need for such an invariant to reconstruct 
sequences of pure states, up to unitary invariance of course. 

Let us first consider the example of a two level system. Let a; be a unit vector 
in R 3 that defines through the Bloch transformation a pure state on the 2x2 
matrices. Let us fix in the canonical basis of C 2 a vector (p x that generates the 
state. Using standard spherical coordinates x = (cos <fi sin 9, sin sin 9, cos 9) 
we choose 

(e^Tilf)' ^(0,tt), 0G(O,2vr(. 

Fix now three unit vectors x, y and z in R 3 . A straightforward computation 
shows that 

. . cos 2 (%0 + cos 2 (^) + cos 2 (%^) - 1 

«.(*(««*)) = 2 2 COB( ^ )cos 2 (¥)COB(¥ 2 ) • 

Here, 9 xy denotes the angular distance between x and y. We need of course 
also sin($(p 3 ,; p y ; p z )) in order to determine <3> uniquely. The angle <£> turns 
out to be zero if and only if there is a permutation (a, (3, 7) of (xy, yz, zx) 
such that 9 a = 9/3 — 7 , i.e. if x, y and z are collinear. 

Consider next the following rather trivial problem. The fidelity between two 
degenerate classical probabilities is either one or zero according to whether 
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the points in which the measures live coincide or not. Suppose that we 
are given for a sequence S = {S 1 ,S 2 , . . . ,S n } of Dirac measures the pairwise 
fidelities Fj k := F(Sj ; 5k). We can obviously reconstruct the sequence 6, up 
to renaming points and reshuffling. It suffices to determine the multiplicity 
of each Sj in the sequence by counting the number of fc's such that Fj k = 1. 
The analogous quantum problem is less trivial. 

Proposition 2. For n = 1, 2, . . . there is generically, up to isometric trans- 
formations and reshuffling, a unique sequence u> = (cui, cu 2 , ■ ■ ■ , ou n ) of pure 
states on a Hilbert space S) with given fidelities Fjk :— F(u)j ; u>k) (j < k) and 
phases $ jk e := ${u)j\ uj k ; uj £ ) (j < k < I). 

Proof. Writing cuj = P^, we shall reconstruct the ipj in a canonical way. 
More precisely 

j 

where {ei, e 2 , . . .} is the standard basis in £ 2 (No), up to multiplication of 
by a complex number zi of modulus one. As tpj is only determined up to a 
phase and using the freedom to choose the phases zi of the standard basis, 
we may assume that Cji > and Cjj > 0. It remains to determine the Cjk 
with 1 < k < j. This can be done recursively. We start with en = I. For 
J = 2 

i i 
C21 = Fi 2 2 and c 22 = (l - F 12 ) 2 . 

i 

For j = 3, c 3 i = F^ 3 . In order to determine the complex coefficient c 32 , we 
use the data F 13 and $ 12 3 

e^ 12a (<p U V 2) (<P2, V?3> (^3, <Pl) = 4 4 4" 

This equation determines (ip 2 , (ps) and from there on c 32 as 

(<P2, V3) = c 2i c 3 i + c 22 c 32 . 
Given the coefficients Ckt with 1 < k < j and 1 < t < k, we determine 
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uniquely the Cje with 1 < I < j by the equations 

91 = 4 

• e~^ 12j (ipi, <p 2 ) {(f2, <Pj) {<Pj, <Pi) 

• e _ ** 13j fa, ip 3 ) (<p 3 , ifj) (ipj, <pi) 




It is already clear from the three state case that only specifying the fidelities 
of pairs in a sequence of pure states is insufficient to reconstruct the sequence. 
A simple parameter count shows that we need (n — l) 2 real parameters for 
a sequence of length n. The fidelities provide n(n — l)/2 parameters and 
we therefore need (n — l)(n — 2)/2 more. The whole set of phases provides 
n(n — l)(n — 2)/6 which is far too much for large n. The argument of the 
proof is however based on the specification of phases with 1 < k < j 
which is the minimal number needed. 

To conclude, we propose an extension of the phase to the mixed state case 
based on the purification of states. To keep things as simple as possible we 
shall only deal with density matrices p on with non-degenerate spectrum 
and trivial kernel. The general case can be dealt with taking the support 
of the density matrix into account. We briefly remind here the purification 
construction, also known as GNS construction. A complex conjugation on $) 
is a complex antilinear transformation ip y-^Tp such that 

Tp = p and (<p, ip) = tp). 

To any operator A on S) we may associate a conjugate operator A by the 
relation Ap> := ATp. It is straightforward to check that A is a complex linear 
operator on S) and that 

aA + B = aA + B, AB = AB : and ~¥ = (A)* . 

Let {/1, / 2 , . . .} be an orthonormal basis of ft diagonalising p: pfj = r,- fj. 
As standard purification of p, we consider the vector 

ttp-=J2 r hj®h m^®i3. (6) 

j 



= F i2 F% F 2 ) 



F 2 



1 1 
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Then p is the marginal of the pure state X i— > (Q p ,XQ p ). The expansion of 
Q p in (6) is the Schmidt decomposition of Q p . There are many pure states 
in the tensor system that have p as a marginal. They are all of the form 
1 <S> U Q p with U a unitary operator on S). We shall use this freedom of 
choice to extend the notion of phase from pure states to mixed ones, namely 
®(pi] P2', Pz) is such that 



e i<S>(pi;P2;ps) _ 1 



inf 

Ui,U 2 ,U 3 



i*(l®£/i Hi ,t®U 2 0,2 ,l®Uz Q 3 ) _ j 



(7) 



Here the infimum is taken over all unitaries on Sj. 

We can rephrase the variational problem (7). Let {gi,g2, . . .} be the orthogo- 
nal basis of eigenvectors of a density matrix a with corresponding eigenvalues 
Sj and let V be another unitary on $), then 

(l<giUQp,l<g)V Qa) — ^2 r j 4 (fj ®U7j,9k® Vg~ k ) 

jk 

= J2 r P^f^9k}(Vg k ,Uf J ) 

jk 

= Tr p*aW*U. 

A study of the phase for mixed states will be the subject of a forthcoming 
paper. 
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